Preview onto version 0.3.4

Here you can do some experiments with Rascal. The computation time is restricted to about 10 secs to prevent high load on server. Have a look at the documentation or at the examples. If you see wrong plots reload may help. You job is cancelled as soon someone else accesses this page, so consider a local installation.

Rascal-toolbox-0.3.4 (C) 2001,2002 Sebastian Ritterbusch

Enter 'help' for help or 'licence' for information about licence and warranty.

 "rascal.rc loaded, predefined e, i, pi"
>"Middle value and higher order interval evaluation";
>
>M(f,Y,z)=f([z,z])+f'(Y)*(Y-z);
>H(f,Y,z)=f([z,z])+f'([z,z])*(Y-z)+f''(Y)*sqr(Y-z)/2;
>
>p=(X-17)^9
 X^9-153*X^8+10404*X^7-412692*X^6+10523646*X^5-178901982*X^4+2027555796*X^3-14772192228*X^2+62781816969*X-118587876497
>
>p(17+1/10)*1.
 1.0000000000E-009
>p(17.1)
   -0.0005035400
>p([17.1,17.1])
 [  -0.0078430176,   0.0070800782]
>
>w=1e-4;Y=[17.1,17.1]+w*[-1,1]/2
 [  17.0999499999,  17.1000500001]
>
>p(Y)
 [-9.1144963283E+007,9.1145229790E+007]
>M(p,Y,mid(Y))
 [-1.2027985001E+003,1.2028008806E+003]
>H(p,Y,mid(Y))
 [  -0.0092670592,   0.0116474509]
>
>y(w)=[17.1,17.1]+w*[-1,1]/2;
>order(t)=[log(diam(p(y(10^t)))) log(diam(M(p,y(10^t),17.1))) log(diam(H(p,y(10^t),17.1)))]/log(10);
>
>order(-4)
 [   8.2607633049    3.3812247688   -1.5371185426]
>
>plotGrid(1);plot(order,-12,0);


>"log-log diagram in width of incoming interval to resulting interval of the three methods, ";
>"whose order can be read from the diagram.";
>

Questions, Problems, Bugs? Please use this form.